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ABSTRACT 

We use the observed abundance and clustering of galaxies from the 2dF Galaxy Red- 
shift Survey to determine the matter density fl m and the linear amplitude of mass 
fluctuations as- We use a method based on the conditional luminosity function, which 
allows straightforward computation of the luminosity dependent bias, 6, of galaxies 
with respect to the matter distribution. This allows us to break the degeneracy between 
bias and erg, which has hampered previous attempts of using large scale structure (LSS) 
data to determine erg. In addition, it allows the inclusion of constraints on the redshift 
space distortion parameter (3 = £l° n 6 /6, and yields average mass-to-light ratios as func- 
tion of halo mass. Using only the luminosity function and the correlation lengths as 
function of luminosity we obtain constraints on fi m and as that are in good agreement 
with COBE. Models with low f2 m and high erg as well as those with high fi TO and low as 
are ruled out because they over (under) predict the amount of clustering, respectively. 
We find the cluster mass-to-light ratio, (M v ; r /L} c i, to be strongly correlated with as- 
Using the additional constraints (M vil /L) c i = (350±70)/i(M/L) Q and (3 = 0.49 ±0.09 
as Gaussian priors significantly tightens the constraints and allows us to break the de- 
generacy between f2 m and as- For flat ACDM cosmologies with scale-invariant power 
spectra we obtain that O ro = 0.271^ and cr 8 = 0.77^.14 (both 95% CL). Adding 
constraints from current CMB data, and extending the analysis to a larger cosmolog- 
ical parameter space, we obtain that fi m = 0.25^}'^ and as — 0.78 ± 0.12 (both 95% 
CL). Thus, we find strong indications that both the matter density £l m and the mass 
variance as are significantly lower than their "standard" concordance values of 0.3 and 
0.9, respectively. We show that cosmologies with il m ~ 0.25 and erg ~ 0.75, as favored 
here, predict dark matter haloes that are significantly less centrally concentrated than 
for the standard ACDM concordance cosmology. We argue that this may solve both 
the problem with the rotation curves of dwarf and low surface brightness galaxies, as 
well as the problem of simultaneously matching the galaxy luminosity function and 
the Tully-Fisher zero-point. 

Key words: cosmology: theory — cosmology: cosmological parameters — galaxies: 
formation — galaxies: halos — large-scale structure of universe — dark matter. 



1 INTRODUCTION 

The most popular cosmological models are variants of the 
cold dark matter (CDM) paradigm within which structure 
grows from adiabatic perturbations imprinted during an 
early inflationary era. These density perturbations therefore 
describe the initial conditions for structure formation, and 
determining the corresponding power spectrum, P(k), is one 
of the holy grails in astrophysics. In this paper we assume 
that P(k) can be well approximated as a power-law (as pre- 
dicted by simple inflationary models), and use the observed 
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abundance and clustering properties of galaxies to constrain 
its normalization. We only consider inflationary CDM cos- 
mologies with adiabatic, scalar-only density perturbations. 
In addition, we assume that neutrinos add a negligible mass 
to the cosmological budget and that the vacuum energy is 
described by a cosmological constant. These cosmologies are 
characterized by 6 parameters: the mass/energy densities 
(in terms of the critical density) of the baryons, f2;,, the cold 
dark matter, tt c and the cosmological constant Oa, the Hub- 
ble parameter h = _ffo/(100kms _1 Mpc -1 ), and the spectral 
index, n s , and normalization, erg, of the initial power spec- 
trum. We define fi m = fit + fl c as the matter density and 
Qk = 1 — fim — as a measure of the spatial curvature. 
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Recent years have seen a tremendous improvement 
in the constraints on these cosmological parameters, with 
a clear concordance cosmology emerging. The location of 
the first peak in the angular power spectrum of cos- 
mic microwave background (CMB) temperature fluctuations 
strongly suggests a flat Universe, for which Qk = (e.g., 
Balbi et al. 2000; Lange et al. 2001; Pryke et al. 2002; Net- 
terfield et al. 2002; Ruhl et al. 2002). An array of different 
observational data, most predominantly that of high redshift 
supernovae (SN) la (e.g., Riess et al. 1998; Perlmutter et 
al. 1999), indicate that Qa ~ 0.7. The contribution from the 
baryons is determined from Big-Bang nucleosynthesis mod- 
els and measurements of the primeval deuterium abundance 
to be Q b h 2 = 0.020 ±0.001 (Buries, Nollett & Turner 2001). 
This is in striking agreement with Ob ft 2 = 0.022 ± 0.004 ob- 
tained directly from recent CMB anisotropy measurements 
(Hanany et al. 2000; de Bernardis et al. 2002; Pryke et 
al. 2002; Netterfield et al. 2002). The HST Key project 
has constrained the Hubble constant to h = 0.72 ± 0.08 
(Freedman et al. 2001). Using this as a prior, the CMB 
data itself yields O m = 0.3 ±0.1 (e.g., Rubifio-Martin et 
al. 2002; Lewis & Bridle 2002), again in excellent agreement 
with the SN la data and with a combined analysis of CMB 
and large scale structure (LSS) data from the 2dF Galaxy 
Redshift Survey (2dFGRS) which yields h = 0.67 ± 0.05 
and Q m = 0.31 ± 0.06 (Percival et al. 2002). Finally, CMB 
anisotropy measurements have constrained the spectral in- 



dex of the initial power spectrum to n s = 1.01_ 



(Pryke et 



al. 2002), in excellent agreement with the inflation paradigm 
which predicts values of n s close to, but not necessarily 
equal to, unity. Shortly after this paper was submitted the 
WMAP-team published their first-year of data on the CMB 
background anisotropies, confirming and strengthening the 
CMB results obtained thus far (Bennett et al. 2003; Spergel 
et al. 2003). 

All these different, mutually consistent, results seem 
to point to a cosmology with (Q m , ^a, h, n 3 ) ~ 
(0.3,0.7,0.04,0.7,1.0), which has become known as the 
"concordance" cosmology. Unfortunately, the normalization 
parameter, as, has not yet been determined with sufficient 
accuracy that it can be included in this list of "concordance" 
parameters. Most attempts to determine erg use either the 
observed abundance of clusters of galaxies, or the cosmic 
shear measured from weak lensing studies. Both methods 
are actually dependent on a combination of both as and 
Q m and virtually all existing constraints are highly degener- 
ate in these two parameters. The cluster abundance method 
has the advantage that it measures the power spectrum di- 
rectly at the scale of interest (i.e., ~ Sh^ 1 Mpc) such that 
no extrapolation is required which depends sensitively on 
the shape of the power-spectrum. The obvious downside of 
this method is that it requires accurate mass estimates of 
individual clusters, which can introduce large uncertainties. 
The weak lensing method has the clear advantage that it 
directly probes the distribution of mass, but with the dis- 
advantage that cosmic shear measurements are extremely 
difficult. Detailed overviews and discussions of the various 
pro and cons of each method, as well as detailed comparisons 
of the results can be found in Jarvis et al. (2003), Smith et 
al. (2002a), Pierpaoli et al. (2002) and references therein. 

For the concordance cosmology (S7 m = 0.3) current es- 
timates for erg range from ~ 0.6 (e.g., Borgani et al. 2001; 



Reiprich & Bohringer 2002; Seljak 2002a; Viana, Nichol & 
Liddle 2002) to ~ 1.0 (e.g., Bacon et al. 2002; Fan & Bahcall 
1998; Pen 1998; Pierpaoli, Scott & White 2001; van Waer- 
beke et al. 2002), an uncertainty that is much larger than 
the typical statistical errors quoted on individual estimates. 
Yet, despite this large uncertainty, the vast majority of nu- 
merical simulations and/or galaxy formation models seem 
to adopt as = 0.9 as the "standard" value. The reason is 
that until fairly recently most studies obtained similar, rea- 
sonably consistent results that are well represented by 



0.5 ±0.05 



(1) 



or a 8 ~ 0.9 for ft m = 0.3 (Edge et al. 1990; Henry & Arnaud 
1991; Bahcall & Cen 1992, 1993; White, Efstathiou & Frenk 
1993; Kitayama & Suto 1996; Eke, Cole & Frenk 1996; Viana 
& Liddle 1996; Pen 1998; Markevitch 1998; Henry 2000). 
Only the more recent estimates listed above have suggested 
significantly different values for as- 

In this paper we present a new method to constrain as 
(and fi m ), based on measurements of the abundance and 
clustering of galaxies. The problem with this kind of LSS 
data is that one needs to transform the observed distribution 
of galaxies to a distribution of mass. This requires knowledge 
of the so-called bias parameter b, conveniently defined as 



<(*) 



P(k) 



(2) 



with -Pgg(fc) the galaxy power spectrum. Unfortunately, the 
bias parameter b depends on scale, galaxy luminosity and 
galaxy type (e.g., Kauffmann, Nusser & Steinmetz 1997; 
Jing, Mo & Borner 1998; van den Bosch, Yang & Mo 2003) 
and is extremely difficult to measure. On sufficiently large, 
linear scales, the scale dependence vanishes (e.g., Pen 1998) 
and the shape of Pgg(fc) is the same as that of the mat- 
ter power-spectrum. Therefore, LSS data is predominantly 
used to constrain the shape parameter fl m h (e.g., Percival 
ct al. 2001; Efstathiou et al. 2002), while constraints on the 
normalization erg require independent measurements of the 
bias b. The method we present here is based on the condi- 
tional luminosity function, introduced by Yang, Mo & van 
den Bosch (2003) and van den Bosch, Yang & Mo (2003), 
and takes this bias (and its luminosity dependence) implic- 
itly into account. Using data from the 2dFGRS only we de- 
rive constraints on Q m and erg that are in excellent agree- 
ment with similar constraints from CMB data. We show 
that additional constraints on the average mass-to-light ra- 
tio of clusters can significantly improve these constraints 
and argue for a flat ACDM cosmology with Q. m ~ 0.25 
and as ~ 0.75. Dark matter haloes in such a cosmology are 
significantly less concentrated than in the standard ACDM 
cosmology with Q m = 0.3 and as ~ 0.9. In fact, we show 
that a small reduction of fi m and/or as, as suggested here, 
significantly alleviates two important problems for the stan- 
dard, as = 0.9, concordance cosmology: the claim that the 
rotation curves of dwarf and low surface brightness (LSB) 
galaxies are inconsistent with CDM haloes, and the failure of 
galaxy formation models to simultaneously match the galaxy 
luminosity function and Tully-Fisher zero-point. 

This paper is organized as follows. In Section 2 we 
present constraints on as and fi m from 2dFGRS data on 
the clustering properties of galaxies. In Section 3 we com- 
bine these data with published data on the CMB and per- 
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Table 1. Cosmological Models. 



ID 


S"2 m 




Jbar 


to 


x 2 (*) 


x 2 (0 


<M vir /L) cl 


P 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A30/90 


0.30 


0.90 


0.14 


13.46 


63.0 


0.91 


512 


0.55 


A30/65 


0.30 


0.65 


0.14 


13.46 


64.6 


1.83 


265 


0.40 


A25/75 


0.25 


0.75 


0.16 


14.16 


63.2 


1.26 


353 


0.43 


A20/65 


0.20 


0.65 


0.20 


15.03 


63.2 


1.29 


268 


0.34 



Parameters of four flat ACDM cosmologies discussed in the text. Column (1) lists the ID, column (2) the matter density, Q m , and 
column (3) the matter power spectrum normalization erg. All four models have Q\ = 1.0 — f2 m , a Hubble constant h = 0.7, a spectral 
index n s = 1.0, and a baryon matter density Q^h = 0.02. The corresponding baryonic mass fraction /t ar = Qt/Q m and age of the 
Universe in Gyrs, to, are listed in Columns (4) and (5), respectively. Columns (6) and (7) list the values of X 2 {<&) (equation [14]) and 
X 2 (?) (equation [15]) of the best fit model for the conditional luminosity function, while the corresponding values for (Mvir/L) c i (in 
h (M/L)q) and f3 (equation [20]) arc listed in columns (8) and (9), respectively. 



form a 6-parameter analysis of flat cosmologies. In Section 4 
we discuss implications of a small reduction in both Q m and 
as with respect to the standard values on the Tully-Fisher 
relation and the rotation curves of LSB galaxies. We sum- 
marize our results in Section 5. 

For the purpose of facilitating the discussion that fol- 
lows, we shall compare four flat ACDM cosmologies with dif- 
ferent f2 m and erg. Table 1 lists a number of characteristics 
of these models, including the baryonic mass fraction and 
the age of the Universe. All four models have Q,i,h 2 = 0.02, 
h = 0.7, and n s = 1.0. Model A30/90 is the standard ACDM 
cosmology with Q m = 0.3 and erg = 0.9. Model A30/65has 
the same concordance value of £7 m but with a lower erg. 
Models A25/75and A20/65 , finally, have both D. m and erg 
lowered with respect to the standard values. 



2 CONSTRAINTS FROM THE ABUNDANCE 
AND CLUSTERING OF GALAXIES 

2.1 The Conditional Luminosity Function 

Yang, Mo & van den Bosch (2003, hereafter YMB03) and 
van den Bosch, Yang & Mo (2003; hereafter BYM03) pre- 
sented a new method to link the distribution of galaxies 
to that of dark matter haloes. This method is based on 
modeling of the conditional luminosity function (hereafter 
CLF), $(L|M)dL, which gives the average number of galax- 
ies with luminosity L±dL/2 that reside in a halo of mass M. 
This CLF is the direct link between the halo mass function 
n(M)AM, specifying the comoving number density of haloes 
of mass M, and the galaxy luminosity function 3>(L)dL, 
specifying the comoving number density of galaxies with lu- 
minosity L, through 



$(L) 



-f 

Jo 



$(L\M)n(M) AM 



(3) 



Alternatively, one may link galaxies to their dark matter 
haloes by specifying the distribution of halo masses associ- 
ated with a given galaxy luminosity (see e.g., Guzik & Seljak 
2002). In fact, both methods are related to each other via 
Bayes' theorem (see eq. [54] in YMB03). 

In CDM cosmologies, more massive haloes are more 
strongly clustered (Cole & Kaiser 1989; Mo & White 
1996, 2002). This means that information on the clustering 



strength of galaxies (of a given luminosity) contains infor- 
mation about the characteristic mass of the halo in which 
they reside. Therefore, an observed luminosity function $(L) 
combined with measurements of the galaxy-galaxy two-point 
correlation function £gg(r, L) as function of luminosity puts 
stringent constraints on $(L|M) (see YMB03 and BYM03). 
In addition, the CLF allows one to compute the average, 
total luminosity of galaxies in a halo of mass M 



(L){M) 



${L\M)LAL 



(4) 



and therewith the average mass-to- light ratios (M j L)(M). 
These can be compared with independent measurements, 
providing further constraints on $(L|M). 

For a given CLF the luminosity function $(L) follows 
directly from equation (3) while, at sufficiently large sepa- 
rations r, the two-point correlation function is given by 



£ gg (r,L) = P(LK dm (r) 



(5) 



Here £dm(0 is the dark matter mass correlation function, 
and b(L) is the average bias of galaxies of luminosity L, 
which derives from the CLF according to 



b(L) 



$(LjM) b(M)n(M) AM. 



(6) 



with b(M) the bias of dark matter haloes of mass M (see 
BYM03 for details). 

The mass function n(M) of dark matter haloes at z — 
can be written in the form 



<M)AM=-^vf{v) 



diner 



dlnM 



AM. 



(7) 



Here p is the mean matter density of the Universe at z = 
0, v = S c /a(M), 5 C is the critical overdensity required for 
collapse at z = 0, f(v) is a function of v to be specified 
below, and o~(M) is the linear rms mass fluctuation on mass 
scale M, which is given by the linear power spectrum of 
density perturbations P(k) as 



\M) 



— f 



P(k) W 2 M {k) k 2 



Ak, 



(8) 
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where Wu(k) is the Fourier transform of the smoothing filter 
on mass scale Aft. 

We use the extended Press- Schechter theory with the el- 
lipsoidal collapse corrections of Sheth, Mo & Tormen (2001) 
and write 



vf(v) = 0.644 



( 1 + ^nr) 



2tt J 



1/2 



exp 



(9) 



and 

b(M) 



1 + 



\fa8 c (z 



•Jaiav 2 ) + y/abiav 2 ) 1 



(av 2 )c + b(l-c)(l-c/2) 



(10) 



with v' = 0.841 v, a = 0.707, b = 0.5 and c = 0.6. The result- 
ing mass function and correlation function of dark matter 
haloes have been shown to be in excellent agreement with 
numerical simulations, as long as halo masses are defined 
as the masses inside a sphere with an average overdensity 
of about 180 (Jing 1998; Sheth & Tormen 1999; Jenkins et 
al. 2001; White 2002). Therefore, in what follows we con- 
sistently use that definition of halo mass when referring to 
M. In Section 2.3 we also define the virial mass, for which 
we use the symbol M v j r . Finally, we use the CDM power 
spectrum P(k) of Efstathiou, Bond & White (1992) with a 
shape parameter 



r = Q, m hcxp(—Q b — V2hflb/Qm) 
(Sugiyama 1995), and compute £dm(r) from 

POO 

£ dm (r) = / A^(fc) 



2 , n sin(fcr) dfc 



kr 



(11) 



(12) 



with A n i(fc) the evolved non-linear power spectrum of the 
dark matter mass distribution, for which we use the fitting 
formula of Smith et al. (2002b). 



2.2 Cosmological Constraints 

Both n(M) and b(M) are cosmology dependent. Whether 
or not a CLF can be found that results in good fits to the 
observed $(L) and £gg( r > L) therefore depends on the cos- 
mological model adopted, and this method may therefore 
be used to place constraints on cosmological parameters. 
In YMB03 we used data from the 2dFGRS to constrain the 
matter density Q, m in a flat ACDM cosmology. We only con- 
sidered cosmological models consistent with the recent weak 
lcnsing constraints of Hoekstra, Yee & Gladders (2002), i.e., 
U8^m 52 = 0.46, and adopted an average cluster mass-to- 
light ratio of (M/L) cl = 500ft (M/L)©. It was shown that 
under these conditions the best-fit model has Qo = 0.3 (and 
thus as — 0.9), in excellent agreement with the standard 
concordance cosmology. 

Here we present a more detailed investigation. Unlike 
in YMB03 we consider erg and Q m as independent model 
parameters and let (M/L) c \ be a free parameter (see Ap- 
pendix A). We start by considering only flat cosmologies 



t Throughout, we adopt a spatial top-hat filter for which 
W M (k) = 3(kR)- 3 [sin(fefl) - fc_Rcos(fc_R)] with M and R related 
according to M = 47rpR 3 /3. 



(Q.K = 0.0) with Q. b h 2 = 0.02, h = 0.7 and n s = 1.0, and 
investigate how the observed &(L) and £gg( r i L) constrain 
Q m and a$. In Section 3 below we relax these assumptions 
and perform a more detailed analysis allowing Qbh , h, and 
n 3 to vary as well. For each (£l m ,ti$) we determine the CLF 
that best fits the 2dFGRS luminosity function of Madgwick 
ct al. (2002) and the amplitude of £ gg (r, L) at the correla- 
tion lengths r (L) obtained from the 2dFGRS by Norberg 
ct al. (2002). Note that throughout this paper all luminosi- 
ties are in the photometric 6j-band (uncorrected for intrin- 
sic absorption by dust), unless specifically stated otherwise. 
We parameterize $(L\M) using a model with 7 free param- 
eters (see Appendix A) and use Powell's multi-dimensional 
direction set method (e.g., Press et al. 1992) to find those 
parameters that minimize 



x 2 = x 2 ($) + x 2 (0 



(13) 



2 /o 

and thus maximize the likelihood Cclf = e x ' . Here the 
first term 



N3, 

E 



(14) 



measures the goodness-of-fit to the observed LF $(ii) with 
errors A$(L»). The second term is defined by 



x 2 (0 = 



i=l 



Z S g(ro(Li), Li) - Z ss (r (Li), Li) 
Ai gg (f (Li),Li) 



(15) 



and measures the goodness-of-fit to the amplitudes of the 
galaxy-galaxy correlation functions. Note that we compute 
£ gg (r, L) at the observed correlation lengths ro(L). These are 
compared to the observed values £ gg (fo(Li),Li), which by 
definition are equal to unity. The errors A£ gg (fo(Li), Li) are 
computed from the errors on fo(L;) and the power-law slope 
of £, gg {r,Li) quoted by Norberg et al. (2002). 

The upper left panel of Figure 1 plots the 68 and 95 
percent confidence levels, obtained by integrating under the 
likelihood surface using the flat priors 0.1 < f2 m < 0.6 and 
0.45 < as < 1.30. As is apparent, the observed abundance 
and clustering properties constrain Q. m and as to a fairly 
broad valley in (f2 m ,ag) parameter space, and all four cos- 
mologies listed in Table 1 are consistent with the obser- 
vations at better than 95 percent confidence. Cosmologies 
with small fl m and large as and those with large O m and 
small as are ruled out at more than 95 percent confidence. 
For comparison, the dashed line indicates the relation be- 
tween as and Q m that best-fits the COBE data'f (Bennett 
et al. 1996). This COBE constraint nicely coincides with the 
valley floor of £cLF(fim, as), therewith indicating that both 
the CMB anisotropies and the galaxy clustering indicate the 
same shape of the power spectrum P(k). 

Understanding how $(!/) and £, gg (r, L) constrain as and 
Q m is somewhat complicated. To guide the discussion, Fig- 
ure 2 plots the bias, b(M), and the evolved, non-linear cor- 
relation function of the dark matter £dm(t*) (equation [12]) 



t For fiat ACDM cosmologies with Q b h 2 = 0.02, h = 0.7, n s = 
1.0 and zero optical depth to rcionization, obtained using the 
parameterization of Liddle et al. (1996) and Bunn & White (1997) 
and the matter transfer function of Eiscnstcin & Hu (1998). 



© 2000 RAS, MNRAS 000, 1-19 



Towards Cosmological Concordance on Galactic Scales 5 





Figure 1. The upper left panel plots the 68 and 95 percent confidence levels of £qlf (fi m , fs) (solid lines). For comparison, the dashed 
line indicates the value of exs as function of £l m obtained from the COBE normalization for a flat ACDM cosmology with fTi^h 2 = 0.02, 
h = 0.7, n s = 1.0 and t = 0.0. The various symbols indicate the four models discussed in the text: A30/90 (solid dot), A30/65 (solid 
square), A25/75 (cross) and A20/65 (solid triangle). The solid and dotted contours in the upper right panel show the 68 and 95 percent 
confidence levels of £cLF(^m, us) when combined with a Gaussian prior on the mass-to-light ratio of clusters (as indicated). The lower 
left panel plots contours of constant (M v ; r /L) c i at values of 250, 350, 450, 1550/i (M/L)q. For clarity, the first three contours are 
labelled. Finally, the lower right panel plots the value of ag that maximizes £cLF(^m, o"s) for given Q m as function of the (M v i r /L) c i of 
the corresponding model. Clearly, higher values of {M v i T /L) c \ imply higher normalizations as, as parameterized by equation (17). The 
light (dark) gray regions indicate the 95 (68) percent confidence levels on the observed (M v i r /L) c i (sec Section 2.3). 



for the four extreme cosmologies considered. First of all, the 
fact that we demand a fit to the observed LF hardly imposes 
any constraints at all; $(£) only depends on the halo mass 
function, and one can always choose an appropriate $(L|M) 
so that one fits perfectly no matter what the shape or 

normalization of n(M). We only include the observed LF as 
a constraint since it sets the normalization of the CLF, which 



allows us to compute mass-to- light ratios. As we show below, 
this proves to be an important asset. The constraints on Q m 
and as are almost solely due to the constraints on £ gg (r, L). 
Typically, increasing as increases the amount of clustering. 
This is immediately apparent from Figure 2 which shows 
that ^DM(r) increases drastically from as = 0.45 to as = 1.3. 
Therefore, in order to keep £gg( r > L) fixed at the observed 
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values the bias b(L) needs to be lowered. This in turn re- 
quires a lower halo bias b(M). However, as is apparent from 
equation (10), b(M) can not become arbitrary small, thus 
imposing an upper limit on as for given f2 m . To emphasize 
the robustness of this result, consider the fl m =0.1 and 
as = 1.3 cosmology. At r = 5ft -1 Mpc, roughly the observed 
correlation length of galaxies with — 18 > Mb J — 51ogft >, — 20 
(Norberg et al. 2002), the dark matter correlation function 
is equal to 4.4. Using that £ 9 <?(r = 5ft -1 Mpc) ~ 1 implies 
that these galaxies require an average bias of b(L) ~ 0.48 
(cf. equation [5]). This, however, is smaller than the mini- 
mum of b(M) (see Figure 2). Therefore, no matter in what 
haloes these galaxies reside, their correlation function will 
always be larger than observed. This cosmology is therefore 
ruled out, a result that is robust against whatever we assume 
regarding the CLF. 

To understand why cosmologies with both high f2 m and 
low as are ruled inconsistent is less trivial. Consider, for 
example, the cosmology with £l m — 0.6 and as = 0.45. 
To explain the galaxy-galaxy correlation lengths obtained 
by Norberg et al. (2002), which range from 4.4ft" 1 Mpc 
to 8.5ft -1 Mpc, requires values for b(L) in the range 1.7 
to 3.2§ . As is apparent from equation 5 and the left panel 
of Figure 2, this can be accomplished by, for example, dis- 
tributing all galaxies over haloes within the relatively narrow 
range lO 13 ft _1 M < M < lO 14 /i _1 M . Therefore, in princi- 
ple, one should be able to find a halo occupation model for 
this cosmology that is perfectly consistent with the data. 
The fact that our model can not accurately match the data 
must therefore reflect a restriction due to the parameteriza- 
tion of the CLF. Although our parameterization has partial 
observational support, results in halo occupation statistics in 
good agreement with semi-analytical models, and is robust 
against small changes (see Appendix B for a detailed dis- 
cussion), we emphasize that at least part of the constraints 
shown in the upper left panel of Figure 1 are due to the par- 
ticular parameterization of the conditional luminosity func- 
tion used. Fortunately, as we shall see below, the lower-left 
corner of the (f2 m , crg)-parameter space is anyway ruled out 
by CMB data, and this restriction therefore does not impact 
on our results. 

2.3 The mass-to-light ratios of clusters 

One of the 7 free parameters of our $(L|M) parameteriza- 
tion is (M/L) c i, the average mass-to-light ratio of clusters of 
galaxies with M > 10 14 ft _1 M o (see Appendix A). Note that 
the halo mass M is defined as the mass within the radius 
inside of which the average overdensity is 180. In order to 
facilitate a more direct comparison with mass-to- light ratios 
available in the literature we convert M/L to M W - 1T /L, where 
Mvir is defined as the mass within the virial radius, inside of 
which the average density is A v i r times the critical density. 
We adopt 

A vir = 18tt 2 + 82(fi m - 1) - 39(fi m - l) 2 (16) 

(Bryan & Norman 1998). To convert M to M v ; r we also need 
a model for the halo concentration as function of halo mass 

§ This is easily verified from equation (5) and the dark matter 
correlation function shown in the right panel of Figure 2. 



and cosmology, for which we use the model of Eke, Navarro 
& Steinmetz (2001). In what follows, whenever masses cor- 
respond to the virial mass we explicitely write M v i r . 

The lower left panel of Figure 1 plots contours of the 
best-fit value of {M v i T /L) c \. Typically, models with higher 
as require higher values of (M v i r /L) c i in order to fit the 
observed abundances and clustering properties of galaxies. 
This is easy to understand; consider once again the Q m = 0.1 
and as = 1.3 cosmology. As discussed in the previous sec- 
tion, this cosmology is ruled out because the dark matter is 
too strongly correlated; there are simply no haloes with suf- 
ficiently small bias, b(M), such that the galaxy-galaxy cor- 
relation function can be made consistent with observations. 
The discrepancy with the data is minimized when all galax- 
ies reside in haloes with M ~ 10 1(, ft _1 M , for which b(M) 
is mimimal. Clearly, galaxies need to avoid more massive 
haloes since otherwise they would overpredict the galaxy- 
galaxy correlation strength even more. This 'avoidance' of 
massive haloes explains the extremely large mass-to-light 
ratio of clusters for this cosmology (see lower left panel of 
Figure 1). If as is lowered, the dark matter becomes less 
strongly clustered. Slowly, galaxies are required to populate 
more and more massive haloes in order to match the ob- 
served correlation lengths. This implies that (M v i v /L) c \ has 
to decrease. In addition, lowering as strongly decreases the 
abundance of massive haloes. Since the abundance of galax- 
ies is fixed by the observed LF, one expects, on average, more 
galaxies per halo. This again will cause M/L to decrease 
with decreasing as- This example shows how constraints on 
the abundance and clustering properties of galaxies yields 
a cluster mass-to-light ratio that depends strongly on both 
fi m and as- 

It is apparent from Figure 1 that (M v i r /L) c i varies 
strongly with position along contours of constant £clf- 
This indicates that we can significantly tighten the con- 
straints on S7 m and as by using additional constraints on 
(M v i r /L) c i. Using a variety of techniques to measure cluster 
masses, Bahcall et al. (2000) obtain (M vir /L B }ci = (330 ± 
77) ft (M/L) q (where Lb is the total B-band luminosity, un- 
corrected for internal extinction). Carlberg et al. (1996), us- 
ing the CNOC Cluster Survey, obtain {M vil /L B ) C \ = (363 ± 
65) ft (M/L) (where we have converted the Gunn r-band 
luminosity to Lb using L r = 1.23Ls, see Bahcall & Comer- 
ford 2002). Taking the average from these two measurements 
we obtain (M vir /L B )ci = (350 ± 70)ft (M/L) . Adding this 
as a (Gaussian) prior to the likelihood £cLF(fim, fs) yields 
the 68 and 95 percent confidence levels shown in the upper 
right panel of Figure 1 (solid lines) . As expected, the prior 
on (M v i r /L) c i drastically tightens the constraints on £7 m and 
as- Marginalizing over as and ft m we obtain fi m = 0. 231^9 
and as = 0.75±°'.i8 

(both 95% CL), respectively. Whereas 
the A30/65, A25/75and A20/65 cosmologies are all consis- 
tent with these constraints at better than 68 percent con- 
fidence, the standard A30/90 cosmology is ruled out at the 
99.99 percent confidence level! In fact, the latter requires a 
cluster mass-to-light ratio of ~ 512ft (M/L) (see Table 1), 
more than 2a away from the observed value. 

In order to quantify the dependence of as on (M v i r /L) c i 
we proceed as follows. From £cLF(£2m, as) shown in the up- 
per left panel of Figure 1, we determine the value of Q m 
that maximizes the likelihood for given as- The lower right 
panel of Figure 1 plots (M v i r /L) c i of this best-fit model as 
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Figure 2. The halo bias b(M) (left panel) and dark matter correlation function £dmM (right panel) for four extreme cosmologies as 
indicated. Note how more massive haloes are more strongly clustered, and that b(M) decreases with increasing erg. This largely counter- 
balances the increase in the clustering strength with increasing erg, thus explaining to a large extent why the observed clustering of 
galaxies only puts very mild constraints on ag (see text for a more detailed discussion). 



function of as (i.e., the relation between as and (M v i T /L) c \ 
along the valley floor) , which is well approximated by 

^^( a^wj ' (17) 

The light (dark) gray areas correspond to the 95 (68) percent 
confidence levels on (M v i r /L) c i. Note that, because contours 
of constant {M v i T /L) c \ are almost independent of Q m (ex- 
cept at Q, m < 0.20), equation (17) is also very similar to the 
maximum as that is consistent with a given (M v ir/L) c i- Re- 
quiring that (M vir /L) c i < 500/i(M/L) Q , for example, trans- 
lates into an upper limit of ag < 0.87. 

Our modeling of the conditional luminosity function 
combined with constraints on the mass-to-light ratio of clus- 
ters has resulted in fairly stringent constraints on Q m . It is 
interesting to compare this with an alternative method to 
constrain Q m from (M v i T /L) c \ directly. One can write 

a m =<^k* cl (is) 

with pl the average luminosity density in the Universe, 
p crit = 3Hi/8nG = 2.78 x 10 n fr 2 M Mpc" 3 the critical 
density, and B c \ the mass-to-light ratio bias of clusters of 
galaxies. Note that this 'bias' is not the same as b defined in 
equation (2). Using (M vil /L) c i = (350 ± 70) ft (M/L) and 
adopting p L = (1.82 ± 0.17) x 10 s h L Mpc~ 3 (Norberg et 
al. 2002), we obtain 

ft m = (0.23 ±0.09) Bci (19) 

which is equal to the estimate of £l m derived above if 
Bci = TO. Remarkably enough this is also the value for 
B c i suggested by observations, which indicate that M/L in- 
creases as a function of scale up to ~ 1/iT 1 Mpc, but then 



flattens out and remains approximately constant (Bahcall, 
Lubin & Dorman 1995; Bahcall et al. 2000). 

2.4 Peculiar velocities 

The galaxy-galaxy correlation function £gg( r ) used above 
is derived from the two-dimensional two-point correlation 
function £(a, 7r), which measures the excess probability over 
random of finding a pair of galaxies with a separation in the 
plane of the sky a and a line-of-sight separation tt. Two pe- 
culiar velocity effects distort this redshift space correlation 
function with respect to the real space correlation function 
£(r) and cause £(a, n) to be anisotropic. On small scales 
the virialized motion of galaxies within dark matter haloes 
cause the correlation function to be extended in the tv di- 
rection (the so called "finger- of- God" effect), while on larger 
scales £(cr, it) is flattened in the n direction due to coherent 
infall of galaxies onto mass concentrations (Kaiser 1987). 
The amplitude of this Kaiser effect depends on the matter 
density Q. m and on the biasing of the galaxy distribution. 
Although the galaxy bias is a complicated function of radial 
scale, galaxy luminosity, and even galaxy type (Kauffmann 
et al. 1997; Jing et al. 1998; YMB03; BYM03), on sufficiently 
large scales the radial dependence should vanish (i.e., the 
bias 6 is linear) and the flattening of £ (a, tt) depends on the 
parameter 

P = SC/b. (20) 

Using data from the 2dFGRS, Hawkins et al. (2002) re- 
cently used a multi-parameter fit to £(cr, n) to obtain f3 — 
0.49 ± 0.09 (68% CL), the most accurate measurement to 
date (cf., Hamilton, Tegmark & Padmanabhan 2000; Taylor 
et al. 2001; Outram, Hoyle & Shanks 2001). As outlined in 
Hawkins et al. (2002), the mean redshift and luminosity of 
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Figure 3. The left panel plots contours of constant (3 computed from the CLF that maximizes Cqlf f° r the given fi m and as- The 
contours shown correspond to (from top to bottom) j3 = 0.67, 0.58, 0.40, 0.31 which are the 68 and 95 percent confidence levels of f} as 
obtained by Hawkins ct al. (2002) from the 2dFGRS. The right panel plots the 68 and 95 percent confidence levels of £cLp(^m,o"g) 
obtained using the priors as indicated. Here G(/3) refers to the Gaussian prior (i = 0.49 ±0.09, while G(/3, T) refers to the combined prior 
of G(/3) plus (M vir /L} c i = (350 ± 70) h (M/L) . 



the galaxies on which this measure is based are z = 0.15 and 
1.4 times the characteristic luminosity L* , respectively. 

The problem with interpreting constraints on /3 in terms 
of constraints on the matter density Q m is that the bias pa- 
rameter b is notoriously hard to measure. However, it can be 
computed straightforwardly from the CLF, and as a func- 
tion of galaxy luminosity. For each of the (fi m , cr 8 )-models we 
compute 6(1. 4L*) from the best-fit CLF using equation (6). 
We convert this bias and £l m to their corresponding values at 
z — 0.15 using the method outlined in BYM03 and compute 
the likelihood Cp that the value of (3 found by Hawkins et 
al. (2002) originates from the given cosmology. The contours 
in the left panel of Figure 3 indicate the corresponding 68 
and 95 percent confidence levels. As is apparent, constraints 
on f3 only yield highly degenerate constraints on Q m and as- 
Note, however, that contours of constant Cp run almost per- 
pendicular to those of constant £clf (cf., upper left panel of 
Figure 1). This indicates that ft and £ g g(r, L), although both 
obtained from £(ct, n), constrain 57 m and as in different ways. 
Using /3 = 0.49 ±0.09 as a Gaussian prior on £cLp(n m , a 8 ), 
results in the 68 and 95 percent confidence levels shown in 
the right panel of Figure 3 (solid contours) . Both fi m and as 
are now well constraint; marginalizing over as and Q. m we 
obtain Q m = 0.27t ;^ and a 8 = 0.87t o jg (both 95% CL), 
respectively. Note that these constraints on Q m and as are 
obtained from the 2dFGRS alone, independent of any other 
data set. 

Hawkins et al. (2002) combined their constraint on 
(3 with a determination of the bias obtained by Verde et 
al. (2002) from an analysis of the 2dFGRS bispectrum. This 
yielded Q. m — 0.23 ± 0.09, consistent with the results pre- 
sented here at the la level. It is extremely reassuring that 
two wildly different methods (CLF modeling versus bispec- 



trum analysis) yield values for the bias in such good agree- 
ment. 

Finally, upon combining the priors on f3 and (M v i r /L) c i 
we obtain the 68 and 95 percent contour levels indicated by 
dotted contours. Marginalizing over as and £l m we obtain 
Q m = 0.27±g;io and a 8 = 0.77±° {° 4 (both 95% CL), respec- 
tively. Whereas the constraint on £l m is perfectly consistent 
with the concordance value, the constraint on the cluster 
mass-to-light ratios strongly favors low-erg models over the 
standard as = 0.9. Keep in mind, however, that this analysis 
is restricted to flat ACDM cosmologies with Q. b h 2 = 0.02, 
h — 0.7, and n s = 1.0. In what follows we relax these con- 
straints and perform a more detailed analysis including ad- 
ditional data on CMB anisotropies. 



3 A JOINT CMB PLUS LSS ANALYSIS 

3.1 Cosmic Microwave Background 

Undoubtedly the most stringent constraints on cosmological 
parameters come from measurements of the angular CMB 
temperature power spectrum. Results from COBE (Bennett 
et al. 1996), BOOMERANG (Netterfield et al. 2002), MAX- 
IMA (Hanany et al. 2000), DASI (Halverson et al. 2002), 
VSA (Scott et al. 2002), CBI (Pearson et al. 2002) and 
WMAP (Bennett et al. 2003) have combined to produce a 
power spectrum of CMB temperature anisotropies with the 
first three acoustic peaks resolved. This has resulted in strin- 
gent constraints on cosmological parameters (e.g., Wang, 
Tegmark & Zaldarriaga 2002; de Bernardis et al. 2002; Pryke 
et al. 2002; Sievers et al. 2002; Efstathiou et al. 2002; Spergel 
et al. 2003) . Unfortunately, because of the unknown amount 
of Thompson scattering experienced by CMB photons, these 
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data only put constraints on the combination a 8 e~ T . Here 
r is the optical depth for Thompson scattering, which de- 
pends, amongst others, on the reionization redshift z la . The 
absence of a clear Gunn-Peterson trough in the spectra of 
high-redshift quasars puts a lower limit on the reionization 
epoch of 2 rc > 5.5 (e.g., Becker et al. 2001; Fan et al. 2001), 
and therewith on r. In addition, the observed height of the 
first acoustic peak requires r < 0.4 (e.g., Griffiths, Barbosa 
& Liddle 1999; Ruhl et al. 2002). Because of these limits on 
t the uncertainty on as is limited to about 30 percent. 

To illustrate the constraints that pre-WMAP CMB data 
put on (f2 m , ^-parameter space we use the Monte-Carlo 
Markov (MCM) chains of Lewis & Bridle (2002; hereafter 
LB02). These consist of large numbers of cosmological mod- 
els, each of which is characterized by N model parameters. 
The chains are constructed using the Metropolis-Hastings 
algorithm, which ensures that the number density distribu- 
tion of these models in the iV-dimensional parameter space 
traces out the posterior probability distribution that the ob- 
served data originates from the given cosmology (see LB02 
and references therein). The clear advantage of this MCM 
chain technique is that it samples the entire iV-dimcnsional 
posterior giving far more information than just the marginal- 
ized distributions of each of the N individual parameters 
(which is what is commonly presented). In addition, it is 
straightforward to extract the marginalized distribution in 
any sub-dimensional parameter space by simple projection. 

Here we use the publicly available N = 6 Markov Chain 
of LB02 consisting of ~ 11500 weakly correlated cosmolog- 
ical models. The six parameters that label each model are: 
flbh 2 , Q c h 2 , h, n s , the redshift z IC at which the reionization 
fraction is a half, and A s , a measure for the amplitude of 
the initial power spectrum. Most of these parameters have 
weak, flat priors (see LB02 for details), except for the Hub- 
ble parameter which has a Gaussian prior of h = 0.72 ±0.08 
(taken from the HST key project, Freedman et al. 2001). All 
models are computed assuming purely adiabatic, scalar-only, 
Gaussian primordial perturbations, Q.k = 0, non-interacting 
CDM, and an effective equation of state for the dark energy 
w = — 1 (corresponding to a pure cosmological constant). 
The data used to construct the chains consists of a combi- 
nation of the results of COBE, BOOMERANG, MAXIMA, 
DASI, VSA, and CBI in the form of band power estimates 
for the temperature CMB power spectrum. The CBI data 
with I > 2000, which is most likely due to non-linear effects, 
is ignored. 

Figure 4 plots the distribution of models of the MCM 
chain in the (Q m , as) parameter space, gray-scale coded ac- 
cording to their value of h (left panel) or r (right panel). As 
mentioned above, the number density distribution of mod- 
els is directly proportional to the marginalized likelihood 
£cMB(n m ,o"8). The two contours indicate the 68 and 95 
percent confidence levels, obtained by integrating under the 
likelihood surface. Marginalizing over as and Q m we obtain 
Q m = 0.26^ i2 and a 8 = 0.80±g;if (95% CL), respectively. 
The 95 percent confidence interval for as therefore covers 
the entire range of values for as quoted in the literature. 
As apparent from the gray-scale coding in the panel on the 
right, this uncertainty on erg is due to the a 8 -e~ T degen- 
eracy mentioned above. Higher values of as (for given Q m ) 
correspond to higher optical depth, and therefore to a higher 
redshift of reionization. The constraints on Q, m are almost 



entirely due to the Gaussian prior on h, as is immediately 
evident from the gray-scale coding in the left panel. Whereas 
the CMB data alone mainly constrains Qk (through the lo- 
cation of the first acoustic peak), the combination of CMB 
data plus a prior on the Hubble constant already puts useful 
constraints on Q. m (e.g., Rubiho-Martin et al. 2002; LB02). 



3.2 Combining CMB plus 2dFGRS 

An extremely attractive feature of the MCM chains is that 
it is straightforward to compute likelihoods including addi- 
tional data. The weight of each model is simply adjusted 
proportional to the likelihood under the new constraint, a 
technique known as importance sampling (see LB02 for de- 
tails). The new, marginalized likelihood is subsequently ob- 
tained from the weighted number density of points as a func- 
tion of flm and as- Using importance sampling we compute 
the combined likelihood £ t ot that both the CMB data plus 
the $(Z/) and £ gg (r, L) originate from the given cosmology. 

The 68 and 95 percent confidence levels of £tot(fim, as) 
are shown in the upper left panel of Figure 5 (solid lines). 
As is apparent from a comparison with the likelihood from 
the CMB data alone (dotted contours), the LSS data adds 
virtually no new constraints. In other words, any model (in 
the 6-dimensional parameter space studied here) that fits 
the CMB data also fits the galaxy correlation lengths as 
function of luminosity. Although the lack of improvement 
on the constraints is perhaps somewhat disappointing, it is 
a beautiful demonstration of the level of concordance among 
two completely independent sets of data. 

As in Section 2 above we can tighten the constraints 
on fi m and as by using priors on the cluster mass-to-light 
ratio and/or (3. Including the Gaussian prior {M v i T /L) c \ = 
(350 ± 70) h (M/L)0 results in the 68 and 95 percent con- 
fidence levels shown in the upper right panel. Clearly, the 
constraints on (M w i T /L) c \ restrict both f2 m and a m to rel- 
atively low values (compared to the parameter space al- 
lowed by the CMB data). Marginalizing over as and Q m 
we obtain Q m = 0.25±g;g| and a 8 = 0.761^2 (both 95% 
CL), respectively. This is inconsistent with the standard 
as = 0.9 at ~ 2a. Fairly similar results are obtained if in- 
stead of the prior on {M v i T /L) c \ we include the Gaussian 
prior P = 0.49 ± 0.09 (fi m = 0.26±g;J| , a 8 = 0.80t°;^, 
see lower left panel of Figure 5). Finally, upon including 
both priors we obtain the confidence levels shown in the 
lower right panel with marginalized probability distributions 
n m = 0.25±g".o? and a 8 = 0.78 ± 0.12 (both 95% CL). 

In summary, the observed clustering properties of galax- 
ies are in perfect agreement with matter power spectra that 
fit the current CMB data. The conditional luminosity func- 
tion models presented here allow us to compute mass-to- 
light ratios and galaxy bias in a completely self-consistent 
way, which in turn allows us to significantly tighten the con- 
straints on £l m and its. The strongest constraints come from 
the observed mass-to-light ratio of clusters, which strongly 
argues for cosmologies with both Q. m and as reduced by 
~ 15 percent with respect to the standard values of 0.3 and 
0.9, respectively. 
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Figure 4. The points correspond to cosmological models from the Monte-Carlo Markov chain of LB02, and are gray-scale coded according 
to their values for h (left panel) and t (right panel). The chain is constructed such that the number density of points is proportional to 
the marginalized likelihood CcMs(^m, &s) with the 68 and 95 percent confidence levels indicates with contours. See text for details. 



4 IMPLICATIONS 

As outlined in Section 1, despite the large uncertainty in 
as, most numerical simulations of structure formation in a 
ACDM cosmology have adopted Q m = 0.3 and as = 0.9. 
This has also been the preferred cosmology for studies of 
galaxy formation. Yet, several studies over the past years, 
including this one, have suggested values for £l m and/or as 
that are reduced with respect to these standard values. 

In this section we investigate the implications that small 
modifications of as and 0, m have on the structure and for- 
mation of galaxies and their associated CDM haloes. In par- 
ticular, we focus on two problems that have been identified 
in recent years for the standard ACDM cosmology; the prob- 
lem of matching the TF zero-point and the inconsistencies 
between observed and predicted rotation curves for dwarf 
and low surface brightness galaxies. 



4.1 The Baryonic Tully-Fisher Relation 

Disk galaxies follow a scaling relation between luminosity 
and rotation velocity known as the Tully-Fisher (TF) rela- 
tion. Since the rotation velocity is a dynamical mass mea- 
sure, the zero-point of the Tully-Fisher relation sets a char- 
acteristic mass-to-light ratio and can therefore, in principle, 
be used to constrain cosmological parameters (see e.g., van 
den Bosch 2000). Here we focus on the so-called "baryonic" 
Tully-Fisher relation between disk mass and rotation veloc- 
ity (McGaugh et al. 2000; Bell & de Jong 2001). Observa- 
tionally, the disk mass is obtained by multiplying the disk 
luminosity with the stellar mass-to-light ratio and adding 



the contribution of the cold gas in the disk. Here we use the 
results of McGaugh et al. (2000)^, who found 

M disk = 1.97xlO^M ( I1 ^- T ) 4 (21) 

Defining dark matter haloes as spheres with an average 
density inside the virial radius, r vir , that is A vir times the 
critical density pcrit, one obtains the following relation be- 
tween the virial mass M V [ T and the circular velocity V v i T at 
the virial radius: 

M vi ^2.33xlO^M ( TU ^) 3 (^^(SS) 

The cosmology dependence enters through A v ; r (eq. [16]). 
The mass of a disk galaxy that forms inside this halo can be 
expressed as 

M disk = / disk M vir (23) 



Here /disk is the fraction of baryonic material in the halo that 
ends up in the disk. Equating (23) to the observed baryonic 
TF relation (21) yields 




This fraction has to obey / d i s k < 1, at least under the stan- 
dard assumption that the baryonic fraction inside dark mat- 
ter haloes can not exceed the universal value. The strongest 



» We have checked that our results are not significantly different 
if we use the somewhat different baryonic TF relation of Bell & 
de Jong (2001). 
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Figure 5. The dotted contours in each of the four panels indicate the 68 and 95 percent confidence levels on Q m and as determined 
from the CMB data plus the Gaussian prior G(h) on the Hubble parameter: h = 0.72 ± 0.08. The solid contours indicate the same 
confidence levels of the combined likelihood £tot with the priors as indicated, using the same nomenclature as in Figure 3. 



constraint on cosmological parameters comes from the high- 
Vwr end. Using (22) and evaluating /disk at M v i r = 3 x 
10 12 /i _1 M Q , which is roughly the maximum halo mass for 
disk galaxies, one obtains 



/ disk = 0.991 hQ 



0.02 



This fraction depends strongly on the ratio of the ob- 
served rotation velocity to the circular velocity of the halo. 
In the CDM paradigm dark matter haloes follow the univer- 
sal NFW density distribution 



p(r) = 



dp 



(r/r s )(l + r/r s ) 



(26) 



(Navarro, Frenk & White 1997). Here r s is a characteristic 
radius, p = f2 m p C rit is the average density of the Universe, 
and S is a dimensionless amplitude which can be expressed 
in terms of the halo concentration parameter c = r v i r /r s as 



5 = 



3 ln(l + c) - c/(l + c) 



(27) 



The circular velocity curve of a NFW density distribution 
reaches a maximum V max at a radius r = 2.163r s . The ratio 
of Umax to the virial velocity is given by 



© 2000 RAS, MNRAS 000, 1-19 



12 van den Bosch, Mo & Yang 




Figure 6. Thin gray linos in both panels indicate contours of constant /disk (equation [25]). The contours corresponding to /disk = 0.1, 
0.2, and 0.3 are labelled for clarity. The gray areas correspond to /disk > 1 (left panel) and C200 > 7.4 (right panel), and indicate the 
regions of parameter space that are excluded according to the observed baryonic Tully-Fisher relation and the observed rotation curves 
of LSB galaxies, respectively (see text for details). The dotted contours in both panels correspond to the 68 and 95 percent confidence 
levels of £oMB(^m^s) (cf. Figures 4 and 5), while the solid contours indicate the same levels of confidence but for Ctot(^m, &s) with 
the Gaussian priors (M v i T /L) c \ = (350 ± 70)h (M/L)q and fi = 0.49 ± 0.09 (cf. lower right panel of Figure 5). The symbols are as in 
Figure 1 and correspond to the four models in Table 1. 



^ = 0.465 , f- ~- . C ..^ ^ (28) 
Kir V ln(l + c) - c/(l + c) v ' 

and is typically larger than unity. We use the model of Eke, 
Navarro & Steinmetz (2001) to compute c as function of halo 
mass and cosmology, and we assume that Kot = Knax, i.e., 
that the observed, flat part of the rotation curve coincides 
with Vmax. Detailed models (Mo, Mao & White 1998; van 
den Bosch 2001; 2002) indicate that in cases were / d i s k is 
sufficiently large V TO t > Knax due to the contribution of the 
disk to the circular velocity. We ignore this effect for the 
moment, but caution that the /disk thus derived may be an 
underestimate. 

In Figure 6 we plot contours of constant /disk as func- 
tion of fi m and as- The gray region in the left panel indi- 
cates the cosmologies for which /disk > 1 and which are thus 
excluded. In addition, we plot the 68 and 95 percent con- 
fidence levels of £cmb (dotted contours) and of £ to t plus 
priors on (M v i r /L) c i and f3 (solid contours). The constraint 
from the TF zero-point (/disk < 1) only rules out models 
which are already inconsistent with the CMB data at more 
than 95 percent confidence. Typically /disk decreases with 
decreasing Q m and decreasing erg. Note that for the stan- 
dard ACDM cosmology, with £l m = 0.3 and as = 0.9, only 
about 30 percent of the available baryons in a 3 x 10 h~ Mq 
halo should be present in the disk. In the case of cosmolo- 
gies favored by the CMB plus LSS data presented here, this 
fraction is even lower, between 10 and 20 percent. Thus, in 
order to explain the zero-point of the baryonic TF relation, 
only a minor fraction of the baryonic mass within r v i r can 
become part of the disk; the remaining baryonic material 
either remains in the halo as hot gas, or is expelled from the 



halo. Note that /disk scales with V v i T such that this fraction 
is even lower in lower mass haloes. This requirement for a 
physical mechanism that can prevent baryons from ending 
up in cold gas or stars is also obvious from a direct compar- 
ison of the halo mass function with the galaxy luminosity 
function (e.g., White & Rees 1978; White & Frenk 1991; 
YMB03) and remains one of the most challenging puzzles in 
the framework of galaxy formation. 

4.2 The Tully-Fisher Zero-Point 

Detailed semi-analytical models for the formation of galax- 
ies have so far been unsuccessful in simultaneously fitting 
the galaxy LF and the TF zero-point. Initially, when these 
studies focussed on Einstein-de Sitter cosmologies, the dis- 
crepancies where found to be extremely large, with mod- 
els tuned to fit the LF predicting much too faint luminosi- 
ties for a given rotation velocity (e.g., Kauffmann, White 
& Guiderdoni 1993; Cole et al. 1994). Heyl et al. (1995) 
showed that this could be significantly improved upon by 
lowering £7 m , but the overall agreement remained unsatis- 
factory. Even for the currently popular ACDM cosmology 
with fi m = 0.3 and as = 0.9 no semi-analytical models 
presented to date has been able to simultaneously fit the 
LF and TF zero-point (Somerville & Primack 1999; Cole et 
al. 2000; Benson et al. 2000, 2002; Mathis et al. 2002). Re- 
markably enough, YMB03 found exactly the same problem 
when comparing the average mass-to-light ratios inferred 
from the CLF with those required to fit the TF zero-point. 
Since the CLF method makes no assumptions about how 
galaxies form, this indicates that the problem is not related 
to the (poorly understood) physics of galaxy formation, but 
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Figure 7. Left panel: The average mass-to-light ratio as function of halo mass for the four cosmological models listed in Table 1, as 
obtained from the best-fit CLF. The solid dot indicates the average mass-to-light ratio of clusters (M v i r /L) c i = 350 ± 70/i (M/L)q , where 
the errorbars indicate the la interval. Right panel: Tully-Fisher relations for the four cosmologies listed in Table 1, obtained from the 
(M v i r /L)(M) shown in the left panel as described in the text. For comparison, the solid dots indicate the Tully-Fisher relation for the 
'local calibrator sample' of Tully & Pierce (2000). 



rather is a problem of more fundamental, cosmological ori- 
gin. 

As shown in Section 4.1, the TF zero-point depends 
strongly on the concentration of dark matter haloes: more 
concentrated haloes have a higher V m ax/V v ir, and thus a 
higher rotation velocity for a given disk luminosity. Since 
halo concentrations are strongly cosmology dependent, the 
TF zero-point problem outlined above may simply indicate 
that the standard cosmological concordance parameters are 
not correct. In order to access the impact of small changes 
in Q m and cr 8 , consider the four cosmologies listed in Ta- 
ble 1. The left panel of Figure 7 plots the average mass- 
to-light ratios, (M v -„/L)(M v - lr ), as obtained from the best- 
fit CLFs. The general trend is the same for all models: 
at M vir < 10 n /i _1 M the mass-to-light ratios strongly in- 
crease with decreasing halo mass. This is required in or- 
der to reconcile the faint slope of the galaxy LF with the 
relatively steep low-mass slope of the halo mass function 
(cf., YMB03). At M vir > 10 14 fe _1 M the mass-to- light ra- 
tio {M/L} is, by construction (see Appendix A) constant 
at (M/L) c \, in good agreement with the observed flatten- 
ing of mass-to-light ratios with increasing scale (Bahcall et 
al. 1995, 2000). Typically, lowering as (while keeping all 
other cosmological parameters fixed) reduces (M/L) c \ and 
increases the mass-to-light ratios on the scales of galaxies 
(cf. models A30/90and A30/65). Simultaneously reducing 
us and fl m along the valley floor of £clf, however, leads 
to a reduction of (M v i r /L) on all mass scales (cf. models 
A30/90 , A25/75 , and A20/65 ). 

We can use the CLFs to compute predictions for the 
TF relation as follows. We assume that TF disk galaxies 
are the brightest galaxies in their haloes. From the CLF, 
it is straightforward to compute the average luminosity of 



the brightest galaxy, {L c }, as function of halo mass (see Ap- 
pendix A), which we convert to magnitudes in the photo- 
metric fej-band using Mq^, = 5.3. Finally, the maximum 
rotation velocity Vmax is obtained from equation (28). As- 
suming that these Knax are equal to the TF rotation veloci- 
ties, we obtain the TF relations shown in the right panel of 
Figure 7. For comparison, we also plot (solid circles) the bj 
band TF relation of the 'local calibrator sample' of Tully & 
Pierce (2000). Here we have converted B-band magnitudes 
to the bj band using bj = B-0.28 (B — V) (Blair & Gilmore 
1982) and adopting B — V — 0.7, which corresponds roughly 
to the average color of disk galaxies (de Jong 1996). 

Compared to the data, the TF relation for the standard 
ACDM cosmology A30/90is somewhat too shallow, clearly 
underpredicting the luminosity of the more massive disk 
galaxies. Since the CLFs are tuned to fit the observed LF, 
this illustrates the TF zero-point problem outlined above. 
Model A30/65 , however, predicts a somewhat steeper TF 
relation, in better agreement with the data, while models 
A25/75 and A20/65 predict TF zero-points that are brighter 
than that of model A30/90by as much as almost an entire 
magnitude. Clearly, relatively small changes in Q. m and/or 
as with respect to the standard concordance values of 0.3 
and 0.9, respectively, strongly alleviate the problem (see also 
Seljak 2002b, c). Although a definite answer requires a more 
thorough investigation, it is extremely encouraging that the 
same cosmological model that is preferred by the observed 
abundance and clustering of galaxies, predicts mass-to-light 
ratios on galactic scales that better match the observed 
TF zero-point. Remarkably enough, as shown in Section 4.1 
above, this cosmology also requires lower values of /disk, and 
thus a stronger efficiency of preventing baryons from ending 
up in the disk. It remains to be seen whether detailed mod- 
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Figure 8. The concentration of dark matter haloes, expressed 
through the dimcnsionlcss quantity Ay/ 2 (equation [31]) as func- 
tion of the maximum circular velocity of the halo. Predictions 
are shown for the four cosmological models listed in Table 1 as 
indicated, and arc computed using the Eke ct al. (2001) model 
for halo concentrations. The gray area labelled 'DATA' indi- 
cates the region occupied by observed galaxies (see Alam et 
al. 2002). Reducing Sl m and/or erg reduces halo concentrations 
and thus A^ 2 (Knas). Note that the Ay^(Vmax) of cosmologies 
A30/65and A25/75 almost exactly overlap. 



els for galaxy formation, can simultaneously match the LF 
and the TF zero-point when considering cosmologies such 
as A25/75. 



4.3 Rotation Curves 

In the past years numerous authors have pointed out that 
the rotation curves of dwarf and low-surface brightness 
(LSB) galaxies indicate dark matter haloes that are less cen- 
trally concentrated than predicted by the standard ACDM 
cosmology with as = 0.9 (e.g., Moore 1994; Burkert 1995; 
van den Bosch et al. 2000; Borriello & Salucci 2001; Blais- 
Ouelette, Amram & Carignan 2001; de Blok, McGaugh & 
Rubin 2001; de Blok & Bosma 2002; Swaters et al. 2003). 

Since the average concentration of a halo of given mass 
is cosmology dependent, a constraint on halo concentrations 
translates into one on cosmological parameters. This prin- 
ciple was recently used by McGaugh, Barker & de Blok 
(2003). Using the high resolution hybrid Ha-HI rotation 
curves of LSB galaxies presented by de Blok, McGaugh & 
Rubin (2001) and de Blok & Bosma (2002) they derived 
an upper limit for the mean observed halo concentration of 
C200 < 7.4. Here C200 = r2oo/r s with r 2 oo the radius inside of 
which the average density of the halo is 200 times the critical 
density. Using the model outlined in Navarro et al. (1997) to 
compute C200 for given cosmology and halo mass McGaugh 
et al. (2003) obtain that 



r . 6 = fi™ /iex P (-n 6 - V2hn b /n m ) - 0.32 (n" 1 - 1) (30) 

This is based on the assumption that all LSB galaxies reside 
in dark matter haloes with M200 = 10 12 /i _1 Mq. Since in 
general one expects that M200 is smaller than this, and since 
c decreases with increasing halo mass, this is a conservative 
upper limit. 

The gray region in the upper right panel of Figure 6 
indicates the part of (S7 m , as) parameter space for which 
haloes with M200 = 10 12 /i -1 M Q have concentrations C200 > 
7.4. Here we have adopted Q,bh 2 — 0.02, h = 0.7, and 
n s = 1.0. As is apparent, the constraint (29) puts stringent 
constraints on Q m and as- In fact, if taken at face value, 
a large fraction of the cosmologies that are consistent with 
the CMB data at the 95 percent confidence level or bet- 
ter is ruled out, including the standard ACDM cosmology 
A30/90 . However, virtually the entire 68 percent confidence 
region obtained from the joined CMB plus LSS analysis pre- 
sented in Section 3 above is consistent with (29). Thus, as 
with the TF zero-point, the problem with the observed ro- 
tation curves of dwarf and LSB galaxies may simply be al- 
leviated by adopting a cosmology with somewhat lower as 
and/or 0, m (see also Zentner & Bullock 2002). 

Unfortunately, the robustness of this result is question- 
able for several reasons. First of all, we have adopted a halo 
mass of M 2 oo = 10 12 /i _1 Mq, whereas in reality LSB galaxies 
will reside in haloes with a variety of halo masses. Secondly, 
in order to obtain a measure of c from an observed rota- 
tion curve one generally fits a mass model to the data. In 
some cases, however, no good fit can be obtained for any 
value of c. This is often interpreted as indicating that dark 
matter haloes do not follow an NFW profile (i.e., the dark 
matter is not cold and collisionless) , but may also indicate 
non-circular motions or other distortions (see e.g., Salucci 
2001; Swaters et al. 2003). 

A more robust comparison of models with data was sug- 
gested by Alam, Bullock & Weinberg (2002). Rather than 
using the halo concentration parameter c, which is diffi- 
cult to extract from an observed rotation curve, Alam et 
al. (2002) introduced the dimensionless quantity 



Ay/2 = 



p{Rv/2) 

pcrit 



1 



v m 



HqR 



V/2 



(31) 



as a more robust measure of the concentration of dark mat- 
ter haloes. Here Rv/2 is defined as the radius at which the 
rotation curve falls to half of its maximum value V max , and 
Av/2 thus measures the mean dark matter density inside 
Rv/2 in units of the critical density p C rit- Unlike c and U v i r , 
the parameters A V / 2 and Umax are easy to obtain from an 
observed rotation curve, without having to resort to mass 
model fitting. Under the assumption that the observed ro- 
tation curve is dominated by the contribution of the dark 
matter, a valid assumption in the case of LSB galaxies, Umax 
and Rv/2 can be converted directly to c and U v i r and thus 
be compared to cosmology dependent predictions. 

For the four cosmologies listed in Table 1 we compute 
Ay/2 as function of Vmax as follows. We assume that CDM 
haloes have a NFW density distribution and that halo con- 
centrations depend on halo mass and cosmology according 
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to the model of Eke et al. (2001). Knax then follows from 
equation (28), while Av/2 is related to the concentration 

parameter c as I 



A v/2 = 3.36 A v 



ln(l + c) - c/(l + c) 



(32) 



Results are shown in Figure 8. The gray area labelled 
'DATA' indicates the region occupied by observed galax- 
ies (see Alam et al. 2002). Models A30/65 , A25/75 , and 
A20/65 predict halo densities inside R V /2 that are a factor 
2, 2, and 4 smaller, respectively, than for the standard con- 
cordance model A30/90 . A reduction in Q. m and/or as with 
respect to the standard concordance values thus strongly al- 
leviates the rotation curve problem. Note that the spread 
in the data is larger than the difference between the four 
cosmologies. However, keep in mind that the model calcu- 
lations reflect the average A v / 2 as function of Knax- Fur- 
thermore, any contribution of the baryons to V max (ignored 
here) may boost Ay^Vmax) of the data with respect to the 
pure dark halo models, and the curves shown here should 
thus be interpreted roughly as lower limits. In addition, we 
point out that the Av/2(Knax) shown here is lower than the 
Av/2(V'max) presented in Alam et al. (2002) for the same 
cosmology, especially at small V max . This is due to the dif- 
ferent models used to compute halo concentrations: whereas 
we rely on the Eke et al. (2001) model, Alam et al. (2002) use 
the model suggested by Bullock et al. (2001) which predicts 
somewhat higher concentrations. These differences indicate 
the theoretical uncertainty regarding the concentrations of 
relatively low mass haloes. Although accurate constraints on 
cosmological parameters are therefore not possible, the gen- 
eral trend is clear: decreasing Q m and/or as reduces the con- 
centrations of dark matter haloes, bringing them in better 
agreement with the data compared to the standard model 
with Q m = 0.3 and as = 0.9. 



5 SUMMARY 

One of the main goals in modern cosmology is to determine 
the initial conditions for structure formation in the early 
Universe, expressed through the initial mass power spec- 
trum P(k). In this paper we assumed that P(k) is a simple 
power-law and used data on the abundance and clustering 
of galaxies to constrain the normalization as- 

Because of the unknown bias of galaxies with respect 
to the mass distribution, previous attempts to constrain 
cosmological parameters from large scale structure (LSS) 
data have mainly focussed on the shape of P(k) rather than 
the normalization. In this paper we analyzed data from the 
2dFGRS using a technique based on the conditional lumi- 
nosity function (introduced by YMB03 and BYM03) that 
self-consistently models the galaxy bias, and its luminosity 
dependence. This method, therefore, allows us to simultane- 
ously constrain both the shape and the normalization of the 
mass power spectrum. In fact, unlike additional methods to 



II This relation between Ay^ and c differs from that in Alam et 
al. 2002 (their cq. 9) and that in Zcntncr & Bullock 2002 (their 
eq. 18), both of which are in error. We are grateful to Andrew 
Zentner for bringing this to our attention. 



constrain as, such as cluster abundances and weak lensing, 
this method allows us to constrain both Q. m and as, rather 
than a combination of both parameters. 

We presented two types of analysis. In the first, we fo- 
cus on flat ACDM cosmologies in which only Q m and as 
are allowed to vary. The other parameters we keep fixed 
at their "concordance" values, i.e., Qbh 2 = 0,02, h — 0.7 
and n s — 1.0. Using the luminosity function 3>(L) and the 
luminosity dependence of the correlation lengths, ro(L), ob- 
tained from the 2dFGRS by Madgwick et al. (2002) and Nor- 
berg et al. (2002), respectively, we obtain constraints on Q. m 
and as that are in excellent agreement with COBE. Mod- 
els with low fl m and high as are robustly ruled out because 
they over-predict the amount of clustering. Models with high 
fi m and low as are also found to be inconsistent with the 
data, but as we argued in Section 2.2, this results mainly 
from our particular parameterization of the CLF. Adding 
the constraint that the quantity f3 = fC 6 /& = 0.49 ± 0.09, 
as obtained by Hawkins et al. (2002) from the redshift space 
distortions in the 2dFGRS, we obtain that Q m = 0.27jlo.i2 
and a 8 = 0.871 o'Je ( both 95% CL )- These constraints, which 
derive only from the 2dFGRS, without any additional data, 
are in good agreement with the constraint 57 m = 0.23 ±0.09 
(68% CL) obtained by Hawkins et al. using constraints on 
the galaxy bias from Verde et al. (2002) based on an analysis 
of the 2dFGRS bispectrum. It is reassuring that such wildly 
different methods yield comparable constraints on cosmo- 
logical parameters and on the bias parameter 6. 

One of the advantages of the CLF models presented 
here is that they allow a straightforward computation of 
the average mass-to-light ratio, {M v i T /L) ci , of clusters of 
galaxies (defined here as systems with masses in excess 
of 10 14 /i _1 Mq). We find a very strong dependence of 
(M v i r /L) c i on a s , which is well parameterized by 



<M vir /i)ci = 350ft (M/L) 



(33) 



Therefore, any additional, independent measurements of the 
mass-to-light ratio of clusters of galaxies allows the con- 
straints given above to be strengthened even further. Tak- 
ing the average value quoted in the literature, {M v i T /L) c i = 
(35O±7O)/i(M/L) (Carlberg et al. 1996; Bahcall et al. 2000) 
we obtain Q m = 0.27^^ and a 8 = 0.771^^4 ( both 95% 
CL). Thus the observed clustering of galaxies, combined 
with constraints on the mass-to-light ratio of clusters, ar- 
gues for a power spectrum normalization that is ~ 2.5a 
lower than the standard value of as = 0.9. This is in agree- 
ment with a rapidly growing number of studies based on 
cluster abundances and cosmic shear measurements. Note 
that our constraint on as mainly owes to the constraint on 
(M v i r /L)ci. Under the assumption that {M v i T /L} c \ is equal 
to the universal mass-to-light ratio, as suggested by the 
fact that M/L is independent of scale on scales larger than 
~ IhT 1 Mpc (Bahcall et al. 1995, 2000), we obtain that 
O m = 0.23 ±0.09 (68 % CL). This is in excellent agreement 
with the constraints on O m given above, thus indicating self- 
consistency. 

In addition to this analysis in which only fi m and as 
were allowed to vary, we also performed a joint analysis of 
the LSS data with pre-WMAP CMB data, this time us- 
ing a 6-parameter analysis of flat ACDM cosmologies. The 
CMB data itself only poorly constrains as because of the 
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well-known degeneracy with the optical depth due to reion- 
ization. Adding the constraints on the galaxy correlation 
lengths does not significantly reduce this degeneracy. In fact, 
any model that fits the CMB data, also fits the observed 
galaxy clustering, strongly suggesting that both data sets 
are consistent with the same matter power spectrum. How- 
ever, including j3 = 0.49 ± 0.09 and (M vir /L) c i = (350 ± 
7O)/i(M/L)0 as Gaussian priors leads to extremely well con- 
strained parameters; Q, m = 0.251q;J° and as = 0.78 ± 0.12 
(both 95% CL). These are consistent at better that the la 
confidence level with the constraints obtained without the 
CMB data using a more restricted set of cosmologies. In 
addition, these results are in perfect agreement with Lahav 
ct al. (2002), who, using a similar combination of 2dFGRS 
and CMB data, derived a 8 = 0.73 ±0.05 (69% CL). Clearly, 
the data argues for a relatively low value of as, and with a 
small preference for a matter density of fi m ~ 0.25 in favor 
of the concordance value of Q. m ~ 0.3. This is in excellent 
agreement with recent work by Melchiorri et al. (2003), who, 
using a combination of CMB and Sloan Digital Sky Survey 
(SDSS) data obtain very similar conclusions. 

Numerous studies have shown that the ACDM concor- 
dance cosmology predicts dark matter haloes that are too 
centrally concentrated. This is apparent from both the ob- 
served rotation curves of dwarf and low surface brightness 
galaxies and from the zero-point of the Tully-Fisher relation. 
However, the majority of these studies have focused on cos- 
mologies with Q m = 0.3 and as = 0.9. Lowering Q m and/or 
as results in less concentrated dark matter haloes. We have 
investigated the effect that small changes in these two cos- 
mological parameters have on the aforementioned problems. 
For a flat ACDM cosmology with Q. m = 0.25 and as = 0.75 
(close to the values preferred by the analysis presented here), 
the halo concentrations are reduced by ~ 20 percent with 
respect to the standard concordance model. This implies av- 
erage densities inside the radius Ry/2, defined as the radius 
where the circular velocity is half the maximum velocity, 
that are a factor 2.5 smaller. Simple tests show that this 
helps significantly in solving both the rotation curve and 
the TF zero-point problem. 
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APPENDIX A: PARAMETERIZATION OF THE 
CONDITIONAL LUMINOSITY FUNCTION 

Following YMB03 and BYM03 we assume that the CLF can 
be described by a Schechter function: 



$(L|M)dL 



t (AY 

L* \L*J 



exp(-L/L*)dL 



(Al) 



Here L* = L*(M), a = a(M) and $* = $*(M); i.e., the 
three parameters that describe the conditional LF depend 
on M. In what follows we do not explicitly write this mass 
dependence, but consider it understood that quantities with 
a tilde are functions of M. 

We adopt the same parameterizations of these three pa- 
rameters as in YMB03, which we repeat here for complete- 
ness. Readers interested in the motivations behind these par- 
ticular choices are referred to YMB03. For the total mass- 
to-light ratio of a halo of mass M we write 



1 ( M \ (M\-^ (M\ 

l) (M) =2 (t)o U) + y 



(A2) 



forM < IO^/i^Mq, while we adopt {M/L)(M) = {M/L) cl 
for M > 10 14 /i _1 Mq. This parameterization has four free 
parameters: two normalizations, (M/L)o and (M/L) c \, a 
characteristic mass Mi, for which the mass-to- light ratio is 
equal to (M/L)o, and one slope 71, which specifies the be- 
havior of (M/ L) at the low mass end. Note that 72 is fixed by 
requiring continuity of (M/L)(M) across M = 10 14 ft _1 M . 

A similar parameterization is used for the characteristic 
luminosity L*: 



M 



L*(M) 
Here 



\Mi) 



(MY 



r(a + 2) 
r(a + i,i)' 



(A3) 



(A4) 



with T(x) the Gamma function and T(o, x) the incomplete 
Gamma function. This parameterization has two additional 
free parameters: a characteristic mass M2 and a power-law 
slope 73. For a(M) we adopt: 
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with Li defined so that a halo of mass M has on average 
one galaxy with L > Li, i.e., 



0.6 0.8 



Figure Al. The dependence of results on small changes in the 
CLF parameterization. Panels on the left and right correspond 
to changes in the parameters M c and Q15, respectively (as in- 
dicated). Upper, middle and lower panels plot x 2 > \°g[{M / L) u] 
and /3, respectively, all as function of the power spectrum normal- 
ization <rg. Sec Appendix B for discussion. 



a(M) = -1.25 + C log(Mi B ) 



(A5) 



Here M15 is the halo mass in units of 10 15 h~ 1 M & , and £ 
describes the change of the faint-end slope a with halo mass. 

Once a and L* are given, the normalization 3>* of the 
CLF is obtained through equation ( A2) , using the fact that 
the total (average) luminosity in a halo of mass M is given 
by 



(L){M) 



f 

Jo 



<fr(L\M) L AL = <3>* L*r(& + 2). 



(A6) 



Finally, we introduce the mass scale M m i n below which the 
CLF is zero; i.e., we assume that no stars form inside haloes 
with M < M min . Motivated by reionization considerations 
(see YMB03 for details) we adopt M min = 10 9 /i _1 M 
throughout. 

This model for $(LjM) thus contains a total of 7 free 
parameters: 2 characteristic masses; M\ and M2, three pa- 
rameters that describe the various mass-dependencies 71, 
73 and and two normalization for the mass-to-light ratio, 
(M/L) and (M/L) cl . 

For the purpose of making predictions for the TF re- 
lation (see Section 4.2), for each halo we define a 'central' 
galaxy whose luminosity we denote by L c . We assume the 
central galaxy to be the brightest one in a halo, consistent 
with the fact that in most (if not all) haloes the brightest 
members reside near the center. The mean luminosity of this 
central galaxy is defined as 



c oc 

J Li 



poo 
J Li 



<t(L\M)dL = 1. 



(A8) 



$(L\M) L AL = i>* L*r(fi + 2,Li/Z*),(A7) 



APPENDIX B: ROBUSTNESS OF RESULTS 

One of the main concerns regarding our constraints on cos- 
mological parameters is the robustness of the results to 
changes in the CLF model. We have introduced two levels of 
parameterization. First of all, it is assumed that <E>(L|M) is 
well fitted by a Schechter form, independent of the halo mass 
M. Secondly, we have assumed various functional forms, 
with a total of 7 free parameters, to describe how the three 
Schechter parameters (L* , a and $*) depend on halo mass. 

We first address the robustness of our results against 
changes in L* , a, and $* by considering modifications in our 
parameterization of (M/L) (equation [A2]). As outlined in 
Appendix A, in the fiducial model we set {M/L) — (M/L) c i 
for haloes with M > M c = 10 14 /i _1 M . Thus it is as- 
sumed that all haloes with M > M c have the same average 
mass-to-light ratio. This is motivated by the fact that var- 
ious studies have suggested that on cluster mass scales the 
M/L varies only weakly with mass (e.g., Bahcall, Lubin & 
Norman 1995; Bahcall et al. 2000; Kochanek et al. 2002). In 
order to investigate the impact of this assumption we com- 
pare results for three different values of M c : 10 14 /i _1 M Q 
(the fiducial value), M c = 10 ls /i _1 M , and M c = 00 (i.e., 
(M/L) oc M 72 for M — > 00). We consider flat ACDM cos- 
mologies with Q, m — 0.3, Qth 2 = 0.02, h — 0.7 and n s — 1.0, 
and compute the best-fit CLFs for a variety of different as. 
The results are shown in the left panels of Fig. Al. The up- 
per panel plots x 2 (equation [13]) as function of as for all 
three models (as indicated). The middle and lower panels 
plot (M/L) 14, the average mass-to-light ratio of haloes with 
M = lO 14 ^ -1 Mq, and (3, both as function of as- Note how 
X 2 (o's) and (3(as) are virtually independent of M c . The only 
quantity that reveals a modest dependency on our assump- 
tion for M c is (M/L) 14 . 

The panels on the right-hand side show the dependence 
of our models to our parameterization of a(M). As is ap- 
parent from equation (A5), in our fiducial model we set 
Qis = q(10 15 /i _1 Mq) = —1.25. This again is motivated 
by observations of the faint-end slope of the LF of clus- 
ters of galaxies (i.e., Beijersbergen et al. 2002; Trentham & 
Hodgkin 2002). The dashed and dotted curves in the right- 
hand panels of Fig. Al correspond to the best-fit CLF mod- 
els with ais = —1.0 and Q15 = —1.50, respectively. Clearly, 
our choise of ais does not have any significant impact on ei- 
ther (M/L)i4 or (3. It does result in small changes of x 2 (°'s) 
but the overal trend remains the same: cosmologies with in- 
termediate values for as are preferred, and cosmologies with 
as Si 1.2 are clearly ruled out. Note also that the absolute 
minimum of \ 2 occurs for ais = —1.25, which is our fiducial, 
observationally motivated, value. 

Together with a number of similar tests described in 
YMB03 and BYM03, Fig. Al indicates that our results are 
robust against (modest) changes in our parameterization of 
L* , a, and 3>*. This leaves, however, the question to what 
extent the assumption of a Schecter function for $(L|M) 
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shapes these results. Our motivation for the Schechter form 
is fivefold: first of all, the (conditional) luminosity function 
of groups and clusters (i.e. systems with M > 10 13 /i _1 Mq) 
is observed to be well fit by a Schechter function (i.e., Tren- 
tham & Hodgkin 2002, Muriel, Valotto & Lambas 1998). 
Second, both the halo mass function, and the (field) galaxy 
luminosity function have the Schecter form, making it the 
natural functional form to choose. Third, since the CLF is 
a reflection of the various physical processes that play a 
role during galaxy formation, the efficiencies of which are 
expected to vary smoothly with halo mass, it seems rea- 
sonable to assume that the CLF, or its functional form, 
does not change abrubtly with halo mass. Fourth, in YMB03 
we presented an alternative form for the CLF and showed 
that this resulted in virtually identical results. Finally, in 
BYM03 we have shown that the halo occupation statistics 
obtained from detailed semi-analytical models of galaxy for- 
mation compare extremely well with those obtained using 
our Schechter-parameterization of the CLF. Nevertheless, 
neither of these arguments conclusively demonstrates that 
our particular parameterization is appropriate over all mass 
scales. Unfortunately, current observational data is not suffi- 
cient to allow completely unparameterized forms of <E>(L|M). 
More work, which we postpone to future papers, is therefore 
required to investigate to what extent alternative functional 
forms for the CLF impact on our results. 
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